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The parton and hadron cascade model PACIAE based on PYTHIA was used to investigate the 
, charged particle elliptic flow in minimum bias pp collisions at the LHC energies. The strings were 

, distributed randomly in the transverse ellipsoid of the pp collision system with major axis of R 

Q\ | (proton radius) and minor axis of R(l — £) before parton rescattering. The charged particle elliptic 

flow as a function of the random number £ and transverse momentum pr were investigated. The 
calculated vije as a function of reaction energy increases monotonously with increasing reaction 
energy up to t/s ~7 TeV and then turns to saturation. With the parton-parton cross section 
enlarges three times in parton rescattering, the rapidity integrated charged particle elliptic flow 



o 



may reach 0.025 at pr ~2 GeV/c in the minimum bias pp collisions at ^J~s~=7 TeV. 
PACS numbers: 25.75.Dw, 24.10.Lx 
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Q\ . Elliptic flow is the 2nd harmonic in fourier expansion of the particle momentum azimuthal distribution [1[ 
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In theoretical studies the elliptic flow is calculated by 
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where O indicates an average of operator O over particles in an event and (O) an average over events [^J . 

The elliptic flow was suggested in hydrodynamic calculations 3] as a signature of collective expansion in the 
relativistic nuclear collisions. The consistency between experimental data of vi(y) as well as V2{pt) at mid-rapidity 
and the corresponding hydrodynamic predictions was regarded as an evidence of the production of partonic matter 
in relativistic nucleus- nucleus collisions ^, The elliptic flow is now one of the most important observables in 
relativistic nucleus-nucleus collisions [6|49|. It has attracted high attention both experimentally and theoretically. 

In non-central collision the reaction zone (overlap region) between two colliding nuclei is spatially asymmetric. 
The rescattering process among the produced particles transfers this spatial asymmetry into the momentum space, 
and the momentum distribution of the produced particles becomes azimuthally anisotropic. This picture is generally 
valid for a large system. In the pp collision at RHIC energy, the average multiplicity (n c h) is n ot big enough for 
collective effects to be detected. However, in the pp collisions at yfs=l TeV, it has already been observed that the 

10 . 111 ]. Therefore the elliptic flow may be measurable in the pp collision at LHC 



multiplicity can reach ^ > 30 
energies and the further theoretical studies are also required. 



The fluctuations in pp interaction region can result in a sizable spatial eccentricity was assumed in jl2|. With 
further assumptions on the nature of this fluctuations and on the eccentricity scaling of elliptic flow, they reported 
that the elliptic flow becomes measurable in high-multiplicity pp collisions at LHC energies. The initial interactions 
produce a number of hot spots was assumed in Their hydrodynamical evolution of two or more hot spots 

generated a sufficiently large elliptic flow. In [y| a Monte Carlo implementation of the dipole (gluon) cascade model 
was introduced to study the elliptic flow in pp collision at LHC energies. They predicted that in J~s= 7 TeV 
pp collision is around 6-7%. The ideal 3+1D hydrodynamic simulation was used to study the V2 in [lq . and they 
concluded that the elliptic flow can occur at least for top multiplicities pp collision at y/s= 14 TeV. In [3| the 
hydrodynamic model calculations showed that if high multiplicity events in the pp collisions at LHC energies contain 
two flux tubes, the elliptic flow may be observable. Many other models also gave their predictions 171 ] . 



We studied the charged particle elliptic flow in minimum bias pp collisions at the LHC energies by a parton 
and hadron cascade model PACIAE 18 1 based on PYTHIA [19] in this letter. It is well known that PYTHIA is 
a model for high energy hadron-hadron (hh) collisions. In the PYTHIA model a hh collision is decomposed into 
the parton-parton collisions. A hard parton-parton collision is described by the lowest leading order perturbative 
QCD (LO-pQCD). The soft parton-parton collision is considered empirically. Because the initial- and final-state 
QCD radiations and multiparton interactions are considered in the parton-parton scattering, the consequence of 
a hh collision is a parton multijet configuration composed of di-quarks (anti-diquarks) , quarks (anti-quarks), and 
gluons, besides a few hadronic remnants. This parton multijet configuration is followed by the string construction 
and fragmentation (hadronization). Therefore one obtains a hadronic final state for a hh (pp) collision. 

For pp collisions the PACIAE model is different from PYTHIA in the addition of the parton rescattering be- 
fore hadronization and the hadron rescattering after hadronization. The PACIAE model consists of the parton 
initialization, parton evolution (rescattering), hadronization, and hadron evolution (rescattering) four stages. 

1. The parton initialization: 

The parton initialization is performed by the PYTHIA model with string fragmentation (hadronization) 
switched-off. One obtains a parton configuration composed of quarks, anti-quarks, and gluons, besides a 
few hadronic remnants for a pp (hh) collision after diquarks (anti-diquarks) being split randomly into quarks 
(anti-quarks). This parton configuration is regarded as quark-gluon matter (QGM) formed in the initial state 
(fireball) of pp collision. 

2. The parton evolution (rescattering): 

The rescattering among partons in QGM is then considered by the 2^2 LO-pQCD parton-parton cross 
sections [gcj]. The differential cross section of a subprocess ij — > kl is 



do; 



ij—^kl 



^-f E . (3) 



dt s . . 

where the factor K (assumed to be 3 in this letter) is introduced considering the higher order and the nonper- 
turbative corrections, a s = 0.47 stands for the effective strong coupling constant, and s, i, as well as u refer to 
the Mandelstam variables. The specific subprocess q\q~2 — > qiq~2, for instance, has 

<3ig2->9i<?2 
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It is regularized by introducing the parton color screen mass /i (=0.63 GeV) as 



V =lf 2+fi2 . (5) 
- _f - 9 (t - p?) 2 



The total cross section of i + j parton collision is then 



fc.i 



— 5 



With the total and differential cross sections above the parton evolution (rescattering) can be simulated by 
the Monte Carlo method until all parton-parton collisions are exhausted (partonic freeze-out). The results 
introduced later were calculated using this pQCD cross sections except special mention. 

3. The hadronization: 

In the hadronization stage, the partonic matter (QGM) formed after parton rescattering is hadronized by the 



Lund string fragmentation regime [19j after st ring reconstruction or by the Monte Carlo coalescence model. The 
former model has been described in detail in [19J. We proposed a phenomenological coalescence model for the 
later one jl8^ . In this model, two partons coalesce a meson and three partons a baryon (antibaryon) according 
to the valence quark structure of hadron and the flavors, positions as well as the momenta of coalescing partons. 
The positions and momenta of coalescing partons are constrained to the phase space requirement 

167T 2 , , h 3 ,_. 

where h 3 /d is the volume occupied by a single hadron in the phase space , d=A refers to the spin and parity 
degeneracies of the hadron, Ar and Ap stand for the position and momentum distances between coalescing 
partons, respectively. In addition, the momenta of coalescing partons have to satisfy the momentum conserva- 
tion. 

4. The hadron evolution (rescattering): 

In this stage the hadronic matter after hadronization proceeds rescattering. It is dealt with by the usual two- 
body elastic and inelastic collisions [21], until the hh collision pairs are exhausted (hadronic freeze-out). The 
rescatterings among n, K, p, n, p(u), A, A, S, 5, 0, J '/y and their antiparticles are considered for the moment. 



The isospin averaged parametrization formula |22|, |23( is assumed for cross section of hh collisions. In addition, 
an assumed constant total cross sections (cr^t = 40 mb, = 25 mb, = 20 mb, and crfJJ = 10 mb) 

and ratio of inelastic to total cross section (0.85) are provided as another option. 



The PYTHIA and PACIAE models with default parameters (except the factor K was assumed to be 3) were used 
to generate pp minimum bias events at y / s=14 TeV. We found that the charged particle V2{pt) from the PYTHIA 
calculations was almost zero and the one from PACIAE was just around 0.0012. Following the idea of hot spots 
[12I [13I lltl ] we assumed that the strings before parton rescattering were distributed randomly in the transverse elliptic 
zone of the pp collision system. The major and minor axes of this ellipsoid were 

a = R, (8) 
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FIG. 1: Charged particle eccentricity and elliptic flow as a function of £ in minimum bias pp collision at y/s =14 TeV. 



and 



RQ--Q, 



(9) 



respectively. In above equations R was the radium of proton and the parameter £ (random number) reflects the 
fluctuation in the initial position distribution in pp collision. The constituent partons in a string were sampled 
randomly in a circle centered at the string and with radius of 1 fm. 
We calculated the charged particle eccentricity 



24] 
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according to parton initial state in position space (later it was simply indicated as e). In this equation a 2 . 



-x 2 



and 



T 2 _ ,,2 



y 2 were the variance of the initial parton distribution in the x and y directions in a given event. 



Figure Q] (a) and (b) give the charged particle eccentricity and the integrated elliptic flow V2 (\f]\ < 2 and < 
Pt < 3GeV/c) as a function of £ in the minimum bias pp collisions at y/s—lA TeV in the PACIAE calculations, 
respectively. Both the eccentricity and elliptic flow increase with increasing £ that may mean the integrated V2 is 
proportional to the initial eccentricity as expected in j^ . 

The charged particle integrated elliptic flow V2 as a function of the transverse momentum px in minimum bias 
pp collisions at -y/i=14 TeV is given in Figj2] In this figure the full circles, squares, and triangles are the results 
calculated with £=0.1, 0.5, and 0.9, respectively. The V2 value increases with increasing £ that may represent the 
final momentum asymmetry is originated from initial position anisotropy. 

Figure [3] shows the reaction energy dependence of scaled elliptic flow «2/e in the minimum bias pp collisions 
calculated with £ = 0.9. In this figure v 2 is integrated over rj and pj>. One sees here that i^/e increases with 
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FIG. 2: Charged particle V2 as a function of pr in minimum bias pp collisions at ^/s=14 TeV. 



increasing reaction energy monotonously up to 7 TeV and then turns to saturation. This behavior is similar to v 2 /e 
varies with centrality observed in relativistic heavy ion collisions reported in 26 1 . 
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FIG. 3: The scaled charged particle elliptic flow, Vz/e, as a function of collision energy in minimum bias pp collisions 



The function of charged particle integrated elliptic flow v 2 versus transverse momentum px varying with parton- 
parton interaction cross section in the minimum bias pp collision at y / s=14 TeV calculated with £=0.9 is shown in 
Fig. 01 The full circles are the results calculated by the pQCD parton-parton cross sections in the parton rescattering 
while the full squares are calculated by three times enlarged parton-parton cross sections. We see in this figure that 
the charged particle integrated elliptic flow increases nearly 60% at ~2 GeV/c if the parton-parton cross section 
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FIG. 4: Charged particle elliptic flow V2 as a function of the transverse momentum pr in minimum bias pp collision at y / s= 
14 TeV with two different parton-parton cross sections. 



is increased by three times. That indicates the importance of the parton rescattering. 

In summary, the parton and hadron cascade model PACIAE was used to study the charged particle elliptic flow in 
the minimum bias pp collisions at LHC energies. The strings were distributed randomly in the transverse ellipsoid of 
the pp collision system with major axis of R (proton radius) and minor axis of i?(l — £) before parton rescattering. 
We investigated the charged particle elliptic flow as a function of the random number £ and transverse momentum 
Pt, respectively. The calculated v^je as a function of reaction energy first increases monotonously with increasing 
reaction energy up to ^fs ^7 TeV and then turns to saturation. We find that the parton rescattering is also very 
important. If the parton-parton cross section enlarges three times, the rapidity integrated charged particle elliptic 
flow V2 increases nearly 60% and reach to 0.025 at px ~2 GeV/c in the minimum bias pp collisions at ^/s=7 TeV. 
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